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Abstract. We investigate the quantum coherent dynamics of an externally driven 
effective two-level system subjected to a slow Ohmic environment characteristic of 
biomolecular protein-solvent reservoirs in photosynthetic light harvesting complexes. 
By means of the numerically exact quasi adiabatic propagator path integral (QUAPI) 
method we are able to include non-Markovian features of the environment and show 
the dependence of the quantum coherence on the characteristic bath cut-off frequency 
LOc as well as on the driving frequency and the field amplitude A. Our calculations 
extend from the weak coupling regime to the incoherent strong coupling regime. In the 
latter case, we find evidence for a resonant behaviour, beyond the expected behaviour, 
when the reorganization energy coincides with the driving frequency. Moreover, we 
investigate how the coherent destruction of tunneling within the two-level system is 
influenced by the non-Markovian environment. 
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1. Introduction 

During tlie past decade, a tremendous advancement in realizing and controlling the 
coherent quantum dynamics in solid-state nanosystems by an external time-dependent 
force has been achieved [1]. The theoretical description of driven quantum systems has 
discerned novel effects related to the control of quantum tunnehng (see, e.g., [2] for a 
comprehensive review) . Classical examples of the achieved degree of quantum control are 
the manipulation of trapped atoms in quantum optics [3] and the control of chemical 
reactions by external laser fields [2]. In quantum optics, it has been experimentally 
demonstrated that the time evolution of a two-level atom can be significantly modified 
by means of a frequency modulated excitation of the atom by use of a microwave field [4j|. 
Moreover, it has been shown that tunneling of an initially localized state in a double well 
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potential can be almost completely suppressed by a properly tailored external driving 
force [5] . In addition, the experimental demonstration of a quantum coherent dynamics 
in a superconducting flux qubit coupled to a superconducting quantum interference 
device (SQUID) has been reported [6]. 

In fact, the coherent dynamics of a quantum system is disturbed by contact with the 
environment |71 [8] . On the one hand, the phase of the quantum coherences is perturbed, 
which leads to decoherence |9], and on the other hand, energy transfer between the 
system of interest and the environment takes place, which gives rise to dissipation [71 [8]. 
Besides the fundamental theoretical interest in such quantum dissipative systems, a 
deep understanding of the various mechanisms leading to dissipation and decoherence 
is of interest in the context of quantum information processing [9[ [10] • In view of the 
requirements expected to be met in the realization of a quantum computer pH], it is of 
the utmost importance to control and stabilize the coherent dynamics of the quantum 
system of interest. Indeed, various candidates for realizing the building blocks of small 
scale quantum information processors with nanoscale solid-state structures have been 
proposed [T2l|13] and also partially realized in ground-breaking experiments [1]. 

To gain a deeper understanding of the quantum dissipative dynamics, a quantitative 
model including the effects of time-dependent driving, decoherence and dissipation is 
required. One generic model to investigate these effects is the time-dependent spin- 
boson model [2[[7j, where the tunnel-sphtting of the coherent two-level system (TLS) is 
usually denoted by A. The environment is commonly described via the spectral density 
J{uj) [7]. In many cases, an Ohmic spectral density, where J{w) oc u, occurs, as in the 
case of an unstructured electromagnetic environment, where all transitions within the 
system are equally damped [?]. 

To take into account the fact that the environmental frequencies are in principle 
limited, a high-frequency cut-off has to be introduced. This frequency scale is 
related to the time scale on which the environmental degrees of freedom evolve [7j. 
In many cases, Uc is chosen to be the largest frequency scale in the problem {ujc ^ A), 
which corresponds to the fact that the environmental fluctuations evolve on the shortest 
possible time scale, and hence the bath is 'fast'. This describes, e.g., electromagnetic 
fluctuations in a crystal host [3 [Ml [15] or in a superconductor [7]. Under this condition, 
a Markov assumption can be made yielding an effective time-local dissipative dynamics. 
This is also the regime in which Bloch-Redfield type master equations are commonly 
used p3], typically in connection with the additional assumption of weak system-bath 
coupling. The opposite case, when c^c <^ A describes an effective adiabatic bath which 
can be treated in a rather simplified manner [7]. There are, however, relevant situations 
when the bath fluctuations occur on a similar time scale on which the system evolves, 
i.e., Uc ~ A. This, for instance, typically occurs in the quantum coherent dynamics at 
the initial stage of photosynthesis in complex biomolecular structures [T6l [T7[ [T8l [T9 l l20] . 
It has only recently been experimentally hinted that the efficiency of the energy transfer 
from the light harvesting antenna complex to the chemical reaction center is promoted 
by the appearance of a quantum coherent dynamics [21], [22]. This hypothesis is 
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further underpinned in pTl [22] , where it has been evidenced that the collective long- 
range electrostatic response of the biomolecular protein environment to the electronic 
excitations should be responsible for the measured long-lived quantum coherences. In 
the experiment reported in [21], a quantum coherent excitonic dynamics in the energy 
transfer among bacteriochlorophyll (BChl) complexes over a time of around 660 fs has 
been measured at a temperature of 77 K. Such a dynamics is highly non-Markovian 
and more elaborate techniques have to be applied in order to provide an appropriate 
theoretical description. 

In what follows, we consider an externally driven TLS subjected to an Ohmic 
environment with an exponential cut-off. Here, the focus is on the non-Markovian 
influence of a finite ujc on the quantum coherent dynamics. The appropriate numerically 
exact tool to investigate non-Markovian dynamics is the quasi-adiabatic propagator 
path integral (QUAPI) method, which we briefly describe in section [2731 We are able to 
investigate the entire parameter regime of weak as well as strong system-bath coupling 
situation beyond the often used scaling limit A <^ Uc- We compare the weak coupling 
approximation based on a Born-Markov approach against the QUAPI method and show 
that the inclusion of non-Markovian effects is indeed necessary to obtain the correct 
result in the regime tUc — A (section l3.ll) . The QUAPI method allows one to go far 
beyond the limiting cases within the Born-Markov approach, since all non-Markovian 
effects are exactly included. 

The investigation of the parameter regime ~ A is motivated by the fact that 
a slow environment is indeed of physical relevance to light-harvesting biomolecular 
complexes which are embedded in a polar solvent [20]. Here, we introduce a model for 
an effective TLS constructed from two interacting chromophores coupled to a protein- 
polar solvent reservoir. Furthermore, we drive the effective biomolecular TLS with an 
external laser and show that the Hamiltonian for the full system can be described in 
terms of the driven spin-boson model. The driven dissipative dynamics is investigated 
with the objective of understanding basic quantum interference phenomena which 
could be realized as proof-of-principle quantum coherent control experiments in light 
harvesting (LH) photosynthetic complexes such as LH II [16] or in artiflcially designed 
TLS nanostructures with speciflc bath properties. Only recently, we have demonstrated 
that in these systems a non-Markovian environment is most successful in generating 
entanglement of two non-interacting qubits which are coupled to the same Ohmic 
environment |23j . 

We compare our results for the undriven TLS with the outcome of real-time Monte 
Carlo simulations [24] and show that the QUAPI approach gives reliable results in 
the non-Markovian strong coupling case (section 13.21) . In addition, we include an 
external time-dependent drive at frequency ui and show in section 13. 3[ for moderate 
and strong driving, that the amplitude of the forced oscillations in the stationary limit 
strongly depends on ui and, moreover, on Uc- Most interestingly, it turns out that a 
slow environment together with a slow drive ui optimizes the forced oscillations in the 
stationary limit. Finally, in section [331 the effect of a slow dissipative environment on 
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the coherent destruction of tunnehng in the TLS is investigated. We find that the bath 
influence is indeed strongest in the scenario c^c — A. 

2. Model and method 

2.1. Model for the dissipative TLS 

The driven two-level system (TLS) bilinearly coupled to a bosonic heat bath is described 
by the generic spin-boson Hamiltonian [7j 

H{t) = Hs{t) + HsB + HB. (1) 

Here, the system Hamiltonian Hs{t) is chosen to be in the basis of the two states |0) and 
|1), each being, for example, the localized charge state of a charge qubit, or the ground 
state and the excited state of a two-level atom. The TLS, with the tunnel splitting A, 
is driven by a time-dependent external driving field of the form e{t) = Acos{uJit) with 
amplitude A and driving-frequency ui yielding 

Hsit) = ^ (Aa, + e{t)a,) , (2) 

with ai=x,z being the Pauli pseudo-spin matrices. The environment to which the TLS is 
bilinearly coupled is modeled as a bath of harmonic oscillators with bosonic creation and 
annihilation operators b\ b and oscillator frequency Ua, hence Hb = J2a ^^abl^ba- The 
coupling between the TLS and the environment is taken into account by the interaction 
Hamiltonian 

HsB = -^(Tz ga {bi + ba) , (3) 

a 

with Qa being the coupling constants. 

2.2. Model for the dissipative photosynthetic light-harvesting effective TLS 

Complex photosynthetic biomolecular structures have recently been shown to exhibit 
quantum interference properties [211 [22] . In particular, energy transfer among the 
excitons within chlorophyll complexes of the sulfur [21j and the purple [22j bacteria 
have provided evidence for long-lived (picosecond time scale) quantum coherent excitonic 
dynamics, a fact that has only recently become associated to the efficiency of the energy 
transfer from the LH antenna complexes to the chemical reaction centres in such large 
biomolecules [211 [221 [23] . 

In this work, we are interested in the dissipative dynamics of the minimal, basic unit 
of a photosynthetic LH complex which would allow us to model and control quantum 
interference mechanisms taking place in such nanostructures. This is done by modeling 
the specific case of an interacting pair of chromophores in a LH II ring [16]. The effective 
single TLS is built up from two chromophores which are coupled by the Forster resonant 
energy transfer A, as sketched in figure [^a), where fvjjj is the transition energy for 
chromophore j . Since the fluorescence lifetime of the single chromophore is much larger 
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than the other time scales of the system [2D], no radiative decay of the excitations is 
taken into account, and we can write the Hamiltonian for the two-chromophore system 
in the basis B2 = {\gi)^\g2) , |5'i)®|e2) , \ei)^\g2) , |ei)(S)|e2)}, where gj (ej) corresponds 
to the ground (excited) state of chromophore j. The A-Forster couphng between the 
two chromophores (figure [TJ^a)) comprises a dipole-dipole interaction which produces a 
non-radiative transfer of an excitation between the chromophores. Such an interaction 
can be written as Hint = ■^('^x'^x + "^y^y)' ^^"^ Hamiltonian of the bare system 
becomes 

Hs = H^® al + al ^ H2 + , (4) 

where is the identity matrix in the space of chromophore i. 

The correlations due to the bath enter through the coupling to a surrounding protein 
environment and to a polar solvent, which, in general, exhibits a frequency dependent 
dielectric constant [25]. For the details of such a mechanism and their possible geometric 
configurations, we refer to [25]. Formally, this process can be modeled by means of 
a quantized reaction field operator i?* = J2a^ai^hci + ^la), where couples the 
chromophores i to the surrounding environment. This, in turn, is modeled as a bath of 
harmonic oscillators which comprise the energy stored in the polar solvent. Such modes 
are represented via the bosonic operators ba, and, as in the previous section, the bath 
Hamiltonian reads Hb = ^„ tvjjJ)^J)a- If ^l^i denotes the change in the dipole moment 
of molecule i during the transition [20], i.e., the difference between the dipole moment 
of the chromophore in the ground and excited states, the two chromophores are coupled 
to their environment via the interaction Hamiltonian 



HsB-- 

The total Hamiltonian for the two-chromophores is then written in the basis B2 as 

H = Hs + Hb + Hsb= (5) 

f -{n+ + v+) \ 

-(^]_ + l^_) 2A 

2A Q^ + V- 

\ n+ + v+ / 

where Q± = ui ± U2, and V± = 6fiiR^ ± (5/i2-R^- Given the biophysical nanostructure 
composition of the LH II rings [H] , we assume that the states of the single chromophores 
couple to the same surrounding protein bath. Consequently we set D = = D'^, and 
drop any subscripts that may differentiate the bath modes associated to chromophores 
1 and 2 in iff 



I A coupling of the two chromophores to two uncorrelated baths can easily be included within the 
QUAPI method, but is not within the aim of the present work. 
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If only singly excited states are taken into account, from the Hamiltonian we 
identify the active environment coupled 2D-subspace {|ei) (g) \g2), \gi) ® |e2)}. In this 
central subspace of ([5]), the effective interacting biomolecular TLS Hamiltonian reads 



(6) 



with ga = Da{Sfii — Sfi2) being the bath coupling constants, and Q = Now the 
defined effective biomolecular TLS has tunneling splitting A; if, as before, such a TLS 
is driven by a time-dependent external driving field e(t) = A cos (uit), and we consider 
that both chromophores have equal transition energies {ui = 002), the Hamiltonian 
becomes equal to the total Hamiltonian H(t) ([1]), and hence we have effectively 
mapped the interacting, environment correlated chromophore system Hamiltonian to 
that of a generic effective spin-boson Hamiltonian. This is esquematically illustrated in 
figure [H^b), where A is the associated "tunneling energy", between the new basis states 
|0) and |1) for the effective biomolecular TLS (formerly the ;B2-basis states \ge) and 
\eg), respectively). 

To gain information on the full dynamics of the system, the initial conditions at 
t = have to be specified. In all our simulations we assume that the full density matrix 
p at the initial time t = factorizes in accordance with 

p{0) = ps{0) ® PB ■ (7) 

Here, ps{0) is the density matrix of the TLS at initial time to = and the decoupled 
bath canonical distribution ps = exp(— //b/Aj^T), where = tr exp(— i/s/fc^^)- 
To keep the bath in thermal equilibrium, the bath is coupled to a not further specified 
super-bath at temperature T. 

The environment is fully characterized by the spectral density J{uj) = '^ada^i^ " 
Ua), being a quasi-continuous function for typical condensed phase applications. It 
determines all bath-correlations that are relevant for the system via the bath auto- 
correlation function [3, [26] 



00 

L{t) = - j dujJ{uj) 



^ fVjJ 

coth cos uot — I sm uot 



(8) 



We emphasize at this point that within our numerical method (see below) the full 
structure of L{t) is included, which is essential to describe non-Markovian features of the 
system dynamics. In particular, we avoid the often used Born-Markov approximation, 
which consists of replacing the real part of L{t) by a 5-function and neglecting its 
imaginary part. 

For what is reported in the following, we use an Ohmic spectral density with an 
exponential cut-off, i.e., 

J{uj) = 2iTauj exp{—u / Uc) , (9) 

where the dimensionless parameter a describes the damping strength and Uc is the cut- 
off frequency. An Ohmic spectral density is a proper choice for, e.g. electron transfer 
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dynamics [HI [27] or biomolecular complexes [121 EDI [231, ^^^^ in the case of 
Josephson junction qubits [28]. In the case of charge qubits subjected to a phonon 
bath, a different spectral density, with a super-Ohmic low-frequency behaviour, results 
better suited [ni[T5]. 

A microscopic derivation for the spectral density associated to the bacteriochloro- 
phylls in the LH II complexes considered in this work has been reported in [25], where 
different forms of a Debye dielectric solvent have been considered. In general, they lead 
to the Ohmic type of spectral density given by ([9]). The dimensionless damping constant 
a of the protein-solvent is directly related to the parameters of the dielectric model [25], 
and has been estimated to be in the range a ~ 0.01 — 1 [20l 125]. The exponential decay 
of the high-frequency cut-off uJc sets the bath characteristic time-scale. If A -C ujc, the 
bath is very fast compared to the effective TLS and loses its memory quickly. Here, a 
Markovian approximation is appropriate and the standard Bloch-Redfield description 
[T6] applies. However, for the considered biomolecular environment, hwc is typically of 
the order of ~ 2 — 8 meV, while the Forster coupling strengths hA ~ 0.2 — 100 meV 
[20} [25]. Hence, the bath responds slower than the dynamics of the excitons evolve and 
non-Markovian effects become dominant, a regime which is accessible only by rather 
advanced techniques. 

Below, we report results in the scaling limit A, T ^ Uc and vary Uc such 
that the system reaches the crossover to the adiabatic limit, i.e., uJc ~ A. Both 
regimes, and the associated crossover, have been studied by real-time Quantum Monte 
Carlo simulations for electron transfer dynamics within the undriven TLS for selected 
parameter combinations [23]. In this work, we go beyond this by including an external 
laser driving and, furthermore, by covering the entire parameter space. Concerning the 
coupling between the TLS and its environment, we study the whole parameter window 
from the weak coupling limit a ^ 1 to the strong coupling limit a ~ 1. 

2.3. The quasi- adiabatic propagator path integral method 

The dynamics of the TLS introduced in the previous section is described in terms of 
the time evolution of the reduced density matrix p{t) which is obtained by tracing over 
the bath degrees of freedom. The TLS dynamics always evolves from the initial state 

P5(0) = |l) (1|. 

In order to investigate the dynamics of the system, we use the QUAPI scheme [29], 
which is a numerically exact iteration scheme which has been successfully adopted to 
many problems of open quantum systems [HI [23], [30] . For details of the iterative scheme 
we refer to previous works [291 EO] and do not reiterate them here. 

The algorithm is based on a symmetric Trotter splitting of the short-time 
propagator K{tk+i,tk) of the full system into a part depending on Hs and Hb + Hsb 
describing the time evolution on a time slice 6t. This is exact in the limit 6t 
but introduces a finite Trotter error to the propagation which has to be eliminated by 
choosing At small enough that convergence is achieved. On the other side, the bath- 
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induced correlations being non-local in time are included in the numerical scheme over a 
finite memory time Tmem = K6t which roughly corresponds to the time range over which 
the bath auto-correlation function L{t) given in ([8]) is significantly different from zero. As 
the environmental fluctuations live on a time scale ~ 1/uic, it is particularly important to 
include the full memory when Uc — A. Note that for any flnite temperature, L{t) decays 
exponentially at long times [7], thus justifying this approach. Moreover, K has to be 
increased, until convergence with respect to the memory time has been found. Typical 
values, for which convergence can be achieved for our spin-boson system, are K < 12 
and 6t ~ (0.1 — 0.2)/A. The two strategies to achieve convergence are countercurrent. 
To solve this, the principal of least dependence has been invoked [30] to flnd an optimal 
time increment in between the two limits K oo and 5t ^ 0. 

3. Dynamics of the driven TLS 

With the time evolution of the reduced density matrix p{t) at hand we can now study 
the dynamics of the driven TLS in terms of the population difference P(t) = {p{t)az) 
of the two states, with the initial condition P{0) = 1. 

3.1. Markovian vs non-Markovian dynamics 

Before addressing the effect of driving, we convince ourselves that the dynamics is indeed 
non-Markovian when Uc — A. For this, we compare the numerical exact QUAPI with a 
Bloch-Redfleld approach. To be speciflc, we compare the QUAPI data with the outcome 
of the weak coupling approximation, see (21.171) and (21.172) in |7j, which results from 
a first order approximation in «. In [31], it has been shown that the outcome of the 
weak coupling approach is equivalent to a Bloch-Redfield treatment. 

In figure [2] (a) the result for a = 0.001 and Uc = lOOA is shown. As expected, the 
agreement between both results is perfect since the system is deep in the Markovian 
(weak coupling) regime. For figure [2] (b) the coupling is increased to a = 0.01 and 
since the bath is still in the scaling limit there is acceptable agreement. In contrast, 
for a cut-off frequency Uc = A, see figure [2] (c), strong deviations between the Bloch- 
Redfield approach and the QUAPI approach arise, indicating that the environment is 
non-Markovian in nature here. Although the smooth cut-off function with the reduced 
spectral weight at frequencies u; ~ c^c > A is fully included within the Bloch-Redfield 
approach, we see strong deviations. This further underpins that, as expected, the Bloch- 
Redfield approach fails when the bath fiuctuates on a time scale comparable to that of 
the system's time evolution. 

When the damping is furthermore increased to a = 0.1 the disagreement between 
the numerical exact QUAPI and the weak coupling approximation is enhanced, see 
figs. [3] (a) and (b). This indicates that the choice of Uc away from the scaling limit 
induces a non-Markovian behaviour of the dissipative TLS dynamics which is further 
underpinned by figure [3] (c), where the relaxation rate jr is shown. We obtain this from 
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a fit of a decaying cosine function witli a single exponential to P{t). In the scaling limit 
A ^ both approaches yield the same 7^, whereas strong deviations occur for smaller 
Uc- Note that ujc enters (21.171) in [7] via the renormalized tunneling amplitude Aeff 
which is included here in the weak-coupling approximation. 

3.2. Undriven case: comparison with Quantum Monte Carlo results 

Next, in order to validate our results in the highly non-Markovian crossover regime 
Uc ~ A, we compare our results with the outcome of real-time Quantum Monte Carlo 
(QMC) simulations |2Yj. In figure H] (main) , the results for the high temperature regime 
are shown, where the parameters are T = and a = 2 (strong coupling). In perfect 
(also quantitative) agreement with [27], P(t) decreases faster for A/uc = 2.4 than 
for A/uJc = 1.6. The cut-off frequency is also related to the reorganization energy 
of the environment [7J, which has the form Ej. = du^^^ = 2ahwc for an Ohmic 
environment. Hence, our results are consistent with the physical expectation, since the 
dynamics of the environment is slower for the smaller ratio A/uc for the same a. For 
the low temperature regime T = O.Auc, shown in the inset of figure HI we also observe 
agreement of the QUAPI results with the outcome of [27 ]. 

Furthermore, we address the case of strong coupling between the TLS and its 
environment. For the scaling limit A <^ cUc it is known that there is a transition at 
a = 0.5 from a coherent decay of P{t) for a < 0.5 to an incoherent decay for a > 0.5 
[3, [32]. To be specific, we choose a = 0.6 at a low temperature T = O.IA, as shown in 
figure [5l Again, this behaviour can be understood in terms of the reorganization energy 
Er, since decreasing Uc here has a similar effect as lowering the damping parameter a. 
By means of our numerical exact method, figure [5] quantifies how the transition between 
coherent and incoherent behaviour depends on Uc for a given a. 

We next put our results in the context of the modeled effective biomolecular TLS. 
For excitations in the LH-II ring of the bacteria chlorophyl molecule (BChls in LH-II 
complexes), the Forster coupling strength ^A ~ 46 — 100 meV, and ~ 2 — 8 meV 
[20[ [25] . and hence the ratio ^ ~ 0.1. For these complexes, a is of the order of 0.1 — 1 
[25] , evidencing a strong coupling between the chromophores and the solvent dielectric. 
If we make a = 0.1, HA = 100 meV, and plot a graph such as the one of figure [S] for 
T = O.lA (~ 116 K), long-lived coherent oscillations are sustained for a time of around 
530 fs {tA ~ 80). Such oscillations can be visualized at fixed Uc = O.IA (not shown). 
This rough estimation is in agreement with the time scale of the coherent oscillations 
recently measured in [2T] for the antenna complex from a green sulphur bacteria that has 
seven BChls per protein subunit. If, on the other hand, we make a = 0.6, uJc = O.IA (as 
shown in figure [5]), coherent oscillations are also found but, due to the stronger coupling 
to the bath, they decay quicker than for the case a = 0.1. We thus observe that for the 
rather simple effective biomolecular TLS model introduced here, which is aimed as a 
guide to the possible realization of further proof-of-principle experiments, we are able to 
demonstrate that the non-Markovian features of the protein-solvent environment help 
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to sustain the quantum coherence mechanisms exhibited by the coupled chromophores 
in a LH-II ring. Since our resuhs are of a general character, in principle derived from 
a generic TLS, we also expect them to be valid in artificially designed nanostructures 
with the specific bath properties described here. 

3.3. Driven case A^O 

We now address the influence of a finite periodic external driving field. The results for 
A = A and ui = 0.05A are shown in figure [61 Similarly to the undriven case, the overall 
decay of P(t) is faster for larger Uc- compare, for example, figure [6] (a) for Uc = 1.5A 
and figure[6] (b) for Uc = 30A, both in the weak-coupling situation a = 0.01. In turn, for 
the case uJc ~ A, the superimposed oscillations due to coherent tunneling survive longer 
than in the scaling-limit A ^ c^c (figure [6](b)), before the system reaches its stationary 
state. There, only the stationary oscillations due to the external driving field survive. 
For a stronger system-bath coupling a the decay of P{t) is faster, as expected, and the 
stationary state is reached faster (not shown). As in the undriven case, the described 
behaviour is qualitatively understandable in terms of the reorganization energy E^-. 

For increasing driving strength to A = lOA, the dependence of the TLS dynamics 
on Uc and a is similar. However, the form of the stationary oscillations turns out to 
be qualitatively different from the case of weak driving A = A. Here, stable stationary 
plateaus emerge, as shown in figured In [1], this has been observed experimentally for 
frequency-modulated excitations of a two-level atom, using a microwave field to drive 
transitions between two Rydberg-Stark states of potassium. In the presence of a slow 
frequency modulation, square wave oscillations of the population difference have been 
detected. They can be understood to mean that the large driving amplitude leads to an 
extreme biasing of the TLS dynamics. The centres of the observed plateaus correspond 
to the extrema of the applied cosine driving field. At the position of these maxima, 
the TLS is maximally biased and since the time-scale of the driving field is much 
smaller than the time-scale of the (unbiased) TLS dynamics due to A, the situation 
of an extreme quasistatic bias results, forming an intermediate self-trapping around the 
maxima of the cosine-like driving field. Indeed, this intermediate self-trapping becomes 
shorter lived, and increasingly washed out, when the cj; ~ A (not shown here). For 
smaller a (figure [7^), the square wave form is superimposed by fast "intrawell" damped 
oscillations which die out on a time scale ~ l/a. For increasing a (figure [7b), the fast 
"intrawell" oscillations are not present anymore and stable plateaus form immediately 
after the "interwell" transition. Due to the strong damping, the level is also strongly 
broadened and the quasistationary picture applies in which the level position is given 
by the mean energy of the broadened level. This is not so much sensitive anymore to 
further variation within a half period of the cosine field. 

In a next step, it is interesting to consider the amplitude A^o of the forced 
oscillations in the stationary limit. In figure [HI we show the results for the nonlinear 
response for the case a = 0.1 and a driving field amplitude A = A. We observe a 
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rather weak dependence on ui when ui < A. Moreover, the response also becomes 
rather weak in the regime of strong detuning. In between these two regimes, we find 
an optimal driving frequency ui ~ A/2, where the amplitude of the forced oscillations 
has a maximum. Note that this rather weak resonance is due to the sizeable damping 
a = 0.1. The maximum can be related to the picture of quantum stochastic resonance 
in an unbiased but strongly damped bistable system ^3\. The response of the system 
becomes optimal when the dissipative tunneling transitions are in synchrony with the 
external drive. In addition to this well known mechanism, we study here the dependence 
on the bath cut-off frequency. The behaviour depicted in figure |8] is, essentially, not 
infiuenced by ujc for this value of the driving strength A and is thus independent of the 
time scale of the environment (up to small quantitative modifications). 

A natural question is whether to expect an enhancement of the response when all 
frequencies are comparable, i.e., tu/ ~ tUc ~ A. This situation is addressed in figure [HI 
where we have chosen ui = 0.5A. We find that for weak to intermediate driving, no 
pronounced resonance appears, as shown in figure [9] (a) for the case A = A. In contrast, 
strong driving can induce a resonant nonlinear response which, however, requires non- 
Markovian dynamics, i.e., Uc ~ A. This is illustrated in figure [9] (b) for A = lOA. A 
slow non-Markovian bath with < A is thus much more efficient in maximizing forced 
oscillations in the stationary limit. This feature occurs for the weak coupling {a = 0.01) 
as well as for the strong coupling case a = 0.6. This constitutes another example where 
a non-Markovian bath helps in increasing coherence in the quantum system [23]. 

In contrast, a Markovian environment in the scaling limit, A -C cUc, largely 
suppresses forced oscillations via its destructive infiuence on coherence. This finding 
is most pronounced in the incoherent strong coupling case, a = 0.6. The dependence 
of Aoo on iUc is again qualitatively understandable via the reorganization energy Ej. in 
the case of weak driving. It explains the reduction of the response for weaker damping 
(a = 0.1 and a = 0.01), compared to the strong damping case a = 0.6 (figure [9] (a)). 
However, for stronger driving, A = lOA, the resonance effect is more pronounced for 
a strong coupling situation, a = 0.6 which cannot be explained in terms of a growing 
reorganization energy Ej.. 

3.4- Coherent destruction of tunneling 

When an isolated symmetric quantum TLS is driven with large frequencies ui ^ A, the 
bare tunneling matrix element effectively becomes renormalized by the zero-th Bessel 
function Jo{x) as A — >• Jo{A/uJi)A = A^s [SlElj. The population of the state |1) with 
an initial preparation |1) = 1 follows as 



This implies that for particular choices of the driving parameters, the Bessel function 
term can be fixed to zero. The first zero then corresponds to A/ui = 2.40482..., and 
then, PL{t) equals unity, since the effective tunnel splitting vanishes. This phenomenon 




(10) 
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is known as coherent destruction of tunneling (CDT); see [2] and references therein for 
further details. 

Naturally, the phenomenon of CDT is influenced when the TLS is coupled to an 
Ohmic environment. A complete standstill of the dynamics will no longer occur, due to 
the relaxation processes induced by the bath. For an Ohmic environment in the scaling 
limit under the assumption of weak damping, this question has been addressed in |35j . 
Here, the QUAPI technique allows us to extend these studies to the case of finite cOc- 
To be specific, we choose a driving frequency ui = 20A. 

For weak coupling, a = 0.01, the CDT is only weakly influenced by the environment, 
as expected, and it turns out that the dependence on the cut-off frequency of the bath 
Uc is also weak, as shown in figure [TO] (left). Nevertheless, we find that the unavoidable 
decay of P(t) occurs the slowest when Uc ~ A. For stronger coupling a = 0.1, as 
shown in figure [10] (right), CDT is more strongly influenced by the magnitude of Uc- 
We observe, again, that a slow bath helps to preserve coherence and the decay of P{t) 
is slow. 

In the regime of strong damping, the situation is different. In fact, we find an 
opposite qualitative behaviour which goes beyond the above given explanation in terms 
of Er- In figure [11], the corresponding results for a = 0.6 are shown. As one can 
see, the decay of P{t) is the fastest when Uc — A. This stems from the fact that for 
uJc ~ A, non-Markovian effects are not negligible, since the time-scales of system and 
environment are of the same order, which is reflected in the fact that the spectral density 
J{uj) has a maximum around A. In this resonant situation, the bath modes around the 
characteristic time scale 1/A of the TLS are most important and disturb the CDT in 
a maximal (resonant) way. Consequently, the decay under the CDT conditions is most 
pronounced. In turn, in the scaling limit, the decay under CDT conditions is slow. 
Interestingly, in both cases of strong coupling {a = 0.1; 0.6) the population difference 
is almost maintained on the considered time windows without decay by the slow bath 
Uc ~ 0.1 A, and the decay of P{t) is the slowest. 

4. Conclusion 

We have investigated the dynamics of the driven spin-boson system in the presence of 
an Ohmic bath. The focus is put on the role of the cut-off frequency Uc- Based on the 
numerical exact QUAPI approach, non-Markovian effects are shown to be relevant when 
UJc ~ A. This effect is more pronounced for strong damping, as expected, and as can be 
seen from the relaxation rate 7^ shown in figure [3] The validity of the QUAPI method 
in the regime ~ A is confirmed by the perfect agreement with published results of 
real time Quantum Monte Carlo simulations [2^ . 

For the unbiased case and for strong couphng, we show that damped coherent 
oscillations exist in the population difference P{t) if the bath has a cut-off frequency 
away from the scaling limit, where the decay is known to be incoherent in nature, as 
shown in figure [5] By comparison with relevant experimental data, we were able to 
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show that our results are directly applicable to biomolecular systems, namely to the 
light harvesting complex LH-II of the bacteria chlorophyl molecule. We have shown 
that a sub-unit of such a biomolecular system can effectively be described by means 
of the spin-boson model and have demonstrated that the non-Markovian features of 
the protein-solvent environment help to sustain the quantum coherence mechanisms 
exhibited within the LH-II complex. Moreover, since our results are based on a model 
of general character, we expect them to apply also for a variety of artificially designed 
nanostructures with the specific bath properties reported here. 

Regarding the coherent control of the effective TLS, we found that a strong 
external driving amplitude in combination with a slow driving laser frequency produce a 
population difference with square wave oscillations in the stationary limit, in agreement 
with the experimental results reported in These square- wave like oscillations stem 
from the fact that the TLS experiences a large quasistatic bias. Moreover, it was shown 
that the amplitude A^o of the forced oscillation in the stationary limit depends strongly 
on the frequency of the driving field. A slow driving cu/ < A field protects the forced 
oscillations against the influence of the dissipative environment, whereas in the case of 
a faster driving, these oscillations are strongly suppressed. This is valid for the weak- 
coupling as well as for the strong coupling case. For larger driving frequencies, the forced 
stationary oscillations are considerably reduced. The stationary amplitude of the forced 
oscillations also strongly depends on the time-scale of the environment determined by 
Uc- For strong driving, we find that the stationary amplitude shows resonant features, 
illustrating that the non-Markovian environment plays a constructive role. 

Finally, we investigated the phenomenon of coherent destruction of tunneling under 
the bath influence. Away from the scaling limit, the influence of the environment is 
weaker and the CDT survives on a significantly longer time scale. For very strong 
coupling, i.e., deep in the incoherent regime, the situation is completely different. Here, 
the CDT is most quickly destroyed when the time scales of the TLS and the environment 
are in resonance. In such a strong coupling scenario, the preservation of coherence 
survives longer in the scaling limit; moreover, the CDT is actually helped by the effect 
of a slow bath (e.g., Uc ~ O.lA), where the decay of P{t) is the slowest. 
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Figure 1. Sketch of the effective fight harvesting biomolecufar TLS formed from a 
pair of A— Forster interacting chromophores. 
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Figure 2. Population difference P{t) for the undriven spin-boson modef A = (sofid 
line: QUAPI, symbols: weak coupling approach), (a) a = 0.001 and cj^ = lOOA 
(Markovian regime), (b) a = 0.01 and iVc = lOOA, and (c) a — 0.01 and ujc — A. The 
temperature is always T = lOA. 
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Figure 3. Same as figure^ but for a strong coupling a = 0.1. (c) shows the relaxation 
rate 7^ as a function of lOc- 
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Figure 4. Population difference P{t) for the TLS for different cut-ofF frequencies cUc 
and a damping parameter a = 2. The temperature T = Aujc (main), and T = OAujc 
(inset). 




Figure 5. Population difference P{t) for the undriven TLS system in dependence on 
the cut-off frequency uic- The temperature T = O.IA, and the damping parameter is 
a = 0.6. 
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Figure 6. Population difference P{t) for the driven TLS. The ampUtude of the driving- 
field is A = A and the driving- frequency is uji = 0.05 A. The temperature T ~ 0.1 A, 
and the damping parameter a = 0.01. The cut-off frequency is (a) Uc = 1.5A and (b) 

ujc = 30A 
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Figure 7. Population difference P{t) for the driven single qubit system. The 
amplitude of the driving-field A = lOA, the driving-frequency loi — 0.05A, the cut-off 
frequency Uc = 0.5A, and the temperature T = O.IA. The damping parameter is (a) 
a = 0.1, and (b) a = 0.6. 
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Figure 8. Amplitude of tlie forced oscillations in the stationary limit Aoo as a function 
of the driving frequency uii, for A = A with a = 0.1 and different cut-off frequencies 
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Figure 9. Amplitude Aoo of the forced stationary oscillations as a function of the cut- 
off frequency lOc for (a) weak to intermediate driving A = A, and (b) strong driving 
A = lOA, for three different values of a. The remaining parameters are uji = 0.5A, 
and T = O.lA. 
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(a) (b) 

Figure 10. Population difference P{t) for the driven TLS under the condition for 
coherent destruction of tunneling, uji — 20A,A = 2.40A. The temperature T — O.lA, 
and the damping is (a) a = 0.01, and (b) a = 0.1. 




Figure 11. Same as figure [TUl but for a strong coupling situation, a — 0.6. 



